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Subterahertz radiation emitted from a variety of short rectangular-, square-, and disk-shaped 
mesas of intrinsic Josephson junctions fabricated from a Bi2Sr2CaCu20s+i single crystal was studied 
from the observed I-V characteristics, far-infrared spectra, and spatial radiation patterns. In all 
cases, the radiation frequency satisfies the conditions both for the ac Josephson effect and for a 
mesa cavity resonance mode. The integer higher harmonics observed in all spectra imply that the 
ac Josephson effect plays the dominant role in the radiation. 
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After the first report of intense continuous tera- 
hertz (THz, ITHz = 10 12 Hz) electromagnetic (EM) 
waves emitted from the intrinsic Josephson junctions 
(IJJs) in the high temperature superconductor (HTSC) 
Bi 2 Sr 2 CaCu208+,5 (Bi-2212) by Ozyuzer et al. [1] with re- 
markably higher intensity than previously generated from 
Josephson junctions [2|, a great deal of interest has been 
drawn not only to the physical mechanism of the radi- 
ation but also to the possible variety of applications in 
the vast fields of science and technology. Nondestructive 
inspections, medical diagnostics, high speed communica- 
tions, imaging technologies for security and defense, etc., 
are potential candidates among them. Presently avail- 
able THz sources such as those using parametric genera- 
tion or pulsed-current methods based on semiconductor 
and/or laser technology are rather weak in output power, 
and are mostly pulsed with incoherent waves. Although 
it is possible to generate output power of even ~ «W in 
a two-dimensional array of Josephson junctions [3j, fre- 
quencies above a few hundred GHz cannot be achieved 
with conventional superconductors due to the small en- 
ergy gap of a few meV (1 meV -£=> 483.5979 GHz). How- 
ever, single crystalline Bi-2212 has an approximately 10 
times larger energy gap, so that in principle it can be em- 
ployed to reach frequencies of several THz. Furthermore, 
it was shown [l| that an entirely new mechanism for THz 
EM wave generation occurs in Bi-2212 due to its highly 
anisotropic layered structure of IJJs |4| . 

Here, we present direct unambiguous evidence that the 
ac Josephson effect is the driving mechanism for the ra- 
diation, in which the fundamental frequency /i equals 
the frequency fj = 2eV/(Nh) of the ac Josephson effect, 
where V is the applied dc voltage across the N IJJs, e 
is the electric charge, and h is Planck's constant. f\ is 
also in resonance with a frequency f mp of a thin EM 
cavity mode, which is inversely proportional to the mini- 
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FIG. 1: (a) Sketch of a cylindrical disk IJJ mesa. The grazing 
angle forming the groove of 10 /im width is approximately 5° . 
(b) Scanning ion microscope images of a rectangular (left), 
square (center), and disk (right) mesa. 



mal mesa cross-section dimension. We demonstrate THz 
(or sub-THz) radiation from mesas of various geometri- 
cal shapes, using rectangular-, square-, and cylindrical 
disk-shaped mesas fabricated by focused ion beam (FIB) 
milling. We first examine the observed radiation frequen- 
cies as functions of the geometrical sizes of the mesas, 
in order to analyze the data by the above two radiation 
conditions [1, 5j. The cylindrical geometry is particularly 
interesting for understanding the radiation, because the 
standing EM cavity mode frequencies /^„ are determined 
by the zeroes of the first derivatives of standard Bessel 
functions, which are incommensurate with the harmonic 
ac Josephson frequency spectrum nfj for integer ro[6(. 
Hence, only one frequency can equal both a cylindrical 
cavity mode frequency and an ac Josephson frequency. In 
contrast, for rectangular mesas, each member of a subset 



TABLE I: List of the parameters of the IJJ mesa samples. See text. 
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FIG. 2: (a) The I-V characteristics (left scale) and output 
radiation intensity detected by the Si-bolometer (right scale) 
from disk mesa D3 at the bath temperature T — 25.0 K of its 
maximal radiation intensity are shown. The inset shows the 
c-axis R-T curve, (b) Details of the shaded high bias region 
in Fig. [21 a) where the emission is observed. 



of higher cavity modes of frequencies /£, can be degen- 
erate with an harmonic of the fundamental ac Josephson 
frequency fj [6J]. For the rectangular mesas used in all 
previous studies [l|, [sj, |7J, it was therefore not possible 
to unambiguously determine the primary mechanism of 
the radiation. Comparisons of the angular distribution 
of the observed radiation with that predicted can also 
provide information as to the relative importance of the 
two mechanisms for the radiation [6j, [8j, |9| . 

The single crystals of Bi-2212 used in the present stud- 
ies were prepared by a floating zone technique [10]. The 
as-grown single crystals were annealed at 650 °C for 24 h 
in argon gas mixed with 0.1% oxygen in order to obtain 
underdoped crystals. The temperature T dependence of 
the c-axis resistance (R-T curve) shows the behavior typ- 
ical of slightly underdoped Bi-2212 crystals with a critical 
temperature T c « 87 K. A small piece of a cleaved crys- 
tal ^50 /im thick was glued onto a sapphire substrate by 
silver paste. Next, silver and gold thin layers were evap- 
orated onto the surface. Then, a groove of width 10 ^m 
and depth of ^2 /j,m was patterned by FIB milling, mak- 
ing an islandlike terrace as sketched in Fig.[Ua). At the 
end of the process, a 10 /im gold wire electrode was con- 
nected to the center of the top metal layer of the mesa by 
silver paste. The sample dimensions and profile curves 
were measured by an atomic force microscope and the 
results are presented in Table |U The cross-sectional pro- 
files are considerably slanted and rounded at the edges, 
resulting in approximately trapezoidal shapes (not shown 



here) [5j, |9(. Typically the top cross-section is approxi- 
mately 10-20% smaller than the bottom one. Thin mesas 
with three different cross-sectional shapes were studied: 
one rectangular, one square, and three circular disks with 
properties listed in Table |TJ The total number N ~ 1000 
of IJ Js is roughly estimated from each mesa height d. 

The I-V characteristics of disk mesa D3, measured by 
sweeping the bias dc current / along the c-axis at the 
bath temperature T = 25.0 K, are displayed in Fig. [2ja) 
together with the radiation intensity detected by the Si- 
bolometer. In Fig. HJb) the radiation region is shown in 
detail in an expanded scale. All mesas that emit THz ra- 
diation studied to date have a critical current density j c 
ranging between 55 and 200 A/cm 2 as listed in Table H] 
The strongest radiation occurs on the return branch of 
the outermost I-V curve. In most cases, the THz radi- 
ation was observed in a very narrow voltage range. For 
example, emission occurs for sample D3 between 0.96 
and 0.98 V, corresponding to / w 10.8 mA. Then, the 
radiation suddenly stops due to a jump to another I- 
V characteristic branch. The radiation power density 
at the detector just before jumping is estimated to be 
1.3 nW/cm 2 , the same order of magnitude as obtained 
previously from rectangular mesas [5| . If the interbranch 
jump did not occur, the radiation intensity would likely 
have grown much stronger [9J]. 

It is also interesting to note that the emission usually 
occurs within a sample-dependent range 6T of the base 
temperature of the mesa, as summarized in Table |T] For 
example, mesa D3 emits between 10 and 50 K. Since the 
constant / f=a 11 mA is fed into mesa D3, it is inevitably 
heated at a rate of about 11 mW, corresponding to the 
enormous heating power density of 8.3 kW/cm 3 . This 
huge heating power density cannot be removed quickly 
enough from the mesa to maintain equilibrium, resulting 
in a considerable rise of the mesa temperature. This local 
heating may induce a chaotic nonequilibrium state and 
may adversely affect the THz radiation Hl,[l2J]. Since the 
heat conduction is progressively worse as T is lowered 13 , 
141 ] and the gap vanishes as T — > T c , these features may 



account for the peculiar temperature dependence of the 
radiation intensity. 

In Fig. El[a), the far- infrared spectra of the THz radia- 
tion from the disk mesas, as measured by a Fourier trans- 
form infrared (FTIR) spectrometer, are shown. In the 
inset to Fig. [21a), the fundamental radiation frequency 



/i = /,/ is plotted as a function of l/(2a), where a is the 
mesa radius. The dashed line represents the resonance 
frequency fj = /ft = xiiCo/(27iVea), where Xh = 1-841 
for the TM (1,1) mode, e is the dielectric constant of Bi- 
2212, and cq is the speed of light in vacuum, predicted by 
the cavity resonance model 15[ , assuming the cylindrical 
cavity boundary condition H^,\ p - a — (dE z /dp)\ p = a = 0. 
The radiation frequency /i is clearly proportional to 
l/(2a). The data were fitted with e = 17.6, in good 
agreement with previous results [a, |9j . Note that this e 
value is about 50% larger than that obtained from in- 
frared spectroscopy [16| . In addition, for each disk mesa, 
the second harmonic at f% = 2/j is clearly visible, and 
these frequencies are easily distinguishable from those 
of the nearest higher disk cavity modes [6J, providing 
unambiguous experimental evidence that the uniform ac 
Josephson current is the primary radiation source. 

In order to further understand the excitation mode in- 
side the disk mesa, the radiation intensity 1 was mea- 
sured at various detection angles 9, relative to the nor- 
mal to the mesa. In Fig. Efb), 1(9) for disk mesa 
D3 is presented. The shadowing effect of the radiation 
from the superconducting Bi-2212 crystal wall outside 
the groove is expected to be negligibly small to first ap- 
proximation. Therefore, the following characteristic fea- 
tures are noted: First, 1(9) is strongly anisotropic, hav- 
ing a maximum around 9 = # max = 20-35° from the 
top (9 = 0°), where a local minimum occurs with inten- 
sity ratio 1(9 = 0°)/I max (6» = 20°) = 0.50-0.65. This 
shallow minimum feature is almost the same as for rect- 
angular mesas, although 6* max is somewhat less than the 
corresponding rectangular mesa value [9|. Second, 1(9) 
rapidly diminishes as 9 approaches 90°. l cav (9) calcu- 
lated [6] by assuming radiation from the TM (1,1) cav- 
ity resonance mode alone is shown by the dashed black 
curve in Fig. [3jb). Clearly, the calculated l cav (9) does 
not fit the experimental data, especially near to 9 = 0°, 
where l cav (9) is a maximum [6|. This disagreement can 
be reduced by introducing a superposition of the radia- 
tion from the uniform ac Josephson current source with 
the same Josephson frequency, as described for rectangu- 
lar mesas [a, |9j . The experimental data are better fitted 
by the dual-source model Xduai(^) with mixing parame- 
ter a = 1.44 [6j, corresponding to 58% of the radiation 
arising from the uniform ac Josephson current source, as 
shown by the solid orange curve in Fig. EJb). 

It is significant that the intensity from the uniform 
source is comparable to that of the fundamental cavity 
mode source, as in rectangular mesas [9[. Since the fun- 
damental cavity mode radiation is enhanced by the cavity 
quality Q-value, a similar enhancement must occur for 
the uniform source radiation 17]. This suggests that the 
radiation from the uniform ac Josephson current in the N 
junctions is coherent, amplifying the output by a factor 
of order N 2 [3| . This interpretation is strongly supported 
by the observation of only integral higher harmonics of 
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FIG. 3: (a) The radiation spectra measured by the FTIR 
spectrometer for three disk mesas. The inset shows the ob- 
served frequencies versus l/(2a). The error bars reflect the 
trapezoidal cross-section profile. The dashed line represents 
the resonance frequency expected from the cylindrical cavity 
TM (1, 1) mode with e = 17.6. (b) Polar plot of the radiation 
intensity for disk mesa D3. Different symbols correspond to 
different runs. The solid orange and dashed black curves are 
the best fit to the dual-source model and its cavity compo- 
nent, respectively. 



the fundamental frequencies shown in Fig. [3ja). Nei- 
ther higher cylindrical cavity excitation frequencies of the 
Bessel type nor subharmonics were observed. The higher 
harmonics are naturally present in the uniform Joseph- 
son current, which radiates coherently, but the nonuni- 
form part of which only can excite one cylindrical cavity 
mode. This experimental evidence clarifies unambigu- 
ously that the THz radiation is mainly generated by the 
uniform mode of the ac Josephson current, in sharp con- 
trast to a number of theoretical predictions [g, Il8l - t20l |. 
Note that although the harmonics arise solely from the 
uniform ac Josephson current source, the larger intensity 
fundamental frequency is a mixture of both the uniform 
ac Josephson and nonuniform cavity sources. Further 
studies of the polarization, coherence, and spatial radia- 
tion patterns of the higher harmonics could provide sup- 
porting information for these conclusions. 

The resonance frequencies of the TM (m,p) modes 
of a rectangular cavity of width w and length L 
are fj = f mp = (c /2jl)^/(m/ L) 2 + (p/wf 1, 
whereas for cylindrical cavities, they are fj — f^ np — 
XmpCo/(27Ty / ea), where Xmp is the pth zero of the deriva- 
tive of J m (z), the standard Bessel function of order m [6J. 



Theoretically, the fundamental modes are the TM (1, 0) 
and TM (1,1) modes for rectangular and cylindrical cavi- 
ties, respectively. Although this is consistent with experi- 
ment for cylindrical mesas, it is very curious that the TM 
(1,0) mode for rectangular mesas was never observed in 
experiment [ll la, l2( • Since the mesa may be considerably 
heated by the dc current, it is very likely that the inho- 
mogeneous heat distribution inside it prevents standing 
EM wave formation, particularly along the longer rect- 
angular dimension. However, such hot spots with T even 
exceeding T c observed by Wang et al. [ll| do not seem 
to be a problem in our experiments, since the fundamen- 
tal frequencies of our disk and square mesas excellently 
obey the linear relation of the cavity resonance frequen- 
cies with 1/a and l/w, respectively. The discrepancies 
in relatively long rectangular mesas with w <§C L — 300- 
400 /urn may have a different origin. In order to check 
this, we fabricated rectangular mesas with smaller L/w 
ratios, as shown in the left and center scanning ion mi- 
croscope images in Fig. [ljb) . The spectroscopic emission 
data from these mesas are presented in Fig. |4j It is ob- 
vious from the inset to Fig. 0] that the frequency at 0.63 
THz of the rectangular mesa Rl with L/w « 3.1-3.3 
obeys the TM (0, 1) cavity resonance mode very well, 
but definitely not the TM (1, 0) resonance mode. These 
rather surprising observations strongly suggest that the 
formation of EM standing waves are restricted by other 
as yet undetermined reasons, and may not be excited 
below some cutoff frequency We propose that this cut- 
off frequency may be the Josephson plasma frequency 
fp = co/(27Ty / eA c ), where A c is the c-axis superconduct- 
ing penetration depth. 
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FIG. 4: The radiation spectra for the square (SI) and rectan- 
gular (Rl) mesas. The inset shows the observed frequencies 
versus the narrower width inverse l/w, where the circles are 
previous data [jj. The black dashed line represents the reso- 
nance frequency expected for the rectangular cavity TM (0, 1) 
mode with e = 17.6. 

In summary, we studied the THz radiation generated 
from IJJ mesas of Bi-2212 with different geometrical 
shapes. Our experimental results clearly demonstrate 
the validity of the cavity resonance model for the fun- 
damental frequency mode of thin square and cylindri- 
cal mesas. The spatial radiation pattern cannot be ex- 



plained by the cavity resonance model alone, but requires 
a substantial contribution from the uniform ac Joseph- 
son current source. More importantly, the frequency 
spectra obtained exhibit higher integral harmonics of the 
fundamental /i, which cannot be obtained from higher 
cavity resonance modes in cylindrical cavities, providing 
unambiguous experimental evidence that the uniform ac 
Josephson current is the primary radiation source. Al- 
though heating effects may significantly alter the I-V 
characteristics, they do not greatly affect the two nec- 
essary radiation conditions: the ac Josephson relation, 
/i = /,/ = 2eV/Nh, and the geometrical cavity reso- 
nance condition fi = /^ or fi = f^, for cylindrical 
or rectangular cavities, respectively, each inversely pro- 
portional to the minimal cross-sectional dimension. We 
propose a further radiation condition that /i > f p , the 
Josephson plasma frequency. 
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